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A generalized Hubbard-Stratonovitch transformation relating an integral over random unitary 
N X N matrices to an integral over Efetov's unitary a-model manifold, is introduced. This trans- 
formation adapts the supersymmetry method to disordered and chaotic systems that are modeled 
not by a Hamiltonian but by their scattering matrix or time-evolution operator. In contrast to 
the standard method, no saddle-point approximation is made, and no massive modes have to be 
eliminated. This first paper on the subject applies the generalized Hubbard-Stratonovitch transfor- 
mation to Dyson's Circular Unitary Ensemble. It is shown how to use a supersymmetric variant of 
' the Harish-Chandra-Itzykson-Zuber formula to compute, in the large- A'^ limit, the n-level correlation 

On ' function for any n. Nontrivial applications to random network models, quantum chaotic maps, and 

lattice gauge theory, are expected. 
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I. INTRODUCTION 

> , 

, The long wave length physics of disordered metals is governed by the so-called diffuson and cooperon modes 

■ describing the coherent long-range propagation of particles and holes. An efficient way to handle the weak-disorder 
^ ^ . perturbation expansion involving these modes is provided by the nonlinear a-model representation pioneered by 

■ Wegner j|] and turned into a nonperturbative tool by Efetov |^ . The great advantage of Efetov's supersymmetric 
\^ ' version of the a model is that it applies not only to diffusive systems but also to systems in the ergodic and quantum 
On ' limits and in the localized regime. In an exciting recent development, the nonlinear a model was extended also to 

[ ballistic systems Q and to deterministic chaotic systems Q. 

As it stands, Efetov's method is applicable to disordered systems that are modeled by a Hamiltonian and where the 
'■^ ' disorder in the Hamiltonian is weak and close to Gaussian. It was not known until very recently how to extend it to 
, a closely related class of models defined by a unitary operator rather than the Hamiltonian. One prominent example 
^ ' are random network models for, say, the random flux problem or the percolation transition in integer quantum Hall 
I systems . These are formulated in terms of random scattering matrices connecting incoming to outgoing channels 
O ' at every node of a network. While there has been some progress in understanding the Hamiltonian (or anisotropic) 
^ ^ limit of these models, the original isotropic models have so far defied analytical treatment. A second example of 
• ^ . importance are periodically driven systems such as the quantum kicked rotor or the quantum kicked top Q . Some 
' of these systems have been found to possess quasi-energy level statistics of a universal type given by Dyson's circular 
random-matrix ensembles. Others exhibit the quantum phenomenon of dynamical localization, which corresponds 
to the Anderson localization of disordered metallic wires, but a rigorous derivation of the correspondence is as yet 
lacking. 

In the present paper I am going to introduce a novel technique, akin to the supersymmetry method of Efetov but 
different from it in its details, which will allow analytical progress to be made with all of the above models. The heart 
of the formalism is a kind of generalized Hubbard-Stratonovitch transformation that takes an integral over random 
unitary N x N matrices (of the type that appears in the Gaussian integral representation of a two-particle Green's 
function for noninteracting electrons) and trades it directly for an integral over Efetov's cr-model space of unitary 
symmetry. The transformation is exact for all N including N = 1, no saddle-point approximation is made, and no 
massive modes need to be eliminated. The present paper is fairly technical as it aims to give a detailed exposition 
of the method along with the necessary proofs. The method is illustrated at the simple example of Dyson's Circular 
Unitary Ensemble (CUE). Other, less trivial applications that have already been worked out and are about to be 
published pM, treat the quantum kicked ro tor a nd the network model of Chalker and Coddington. 



The paper is organized as follows. In Sec. H A the generalized Hubbard-Stratonovitch transformation is introduced 
in the context of the calculation of the generating function for spectral correlations of the CUE. An integral over 
U(iV) is replaced by an integral over a pair of supermatrices Z, Z . The geometric meaning of these supermatrices 



is elucidated in Sec. II B, where it is argued that they are coset representatives for Efetov's coset space G/H of 



unitary symmetry. In the Appendix it is shown that the matrix elements of Z can be viewed as the parameters of a 
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generalized coherent state obtained by acting with Efetov's symmetry group G (or, rather, its complexification) on 
a certain lowest-weight module. By embedding this module into a Fock space of bosons and fermions we succeed in 
giving an easy proof of the generalized Hubbard-Stratonovitch transformation. In Sec. II C the n-level correlation 
function of the CUE is calculated for all n, by using a supersymmetric version of the Harish-Chandra-Itzykso n-Zu ber 
formula. Dyson's Circular Orthogonal Ensemble and the Circular Ensemble of type C are discussed in Sees. 
IV , demonstrating the versatility of the new method. A summary and outlook is given in Sec. 0. 

The present work was inspired by unpublished results of Alexander Altland, who had developed a similar (but less 
useful) Hubbard-Stratonovitch scheme. I am grateful to him for many discussions of this subject. 



Ill 



and 



II. CIRCULAR UNITARY ENSEMBLE (TYPE A) 

Consider U(-/V), the unitary group in N dimensions. This group, as any unimodular Lie group, comes with a natural 
integration measure, the so-called Haar measure dU, which is invariant under left and right translations U t-^ UlUUr. 
We normahze dU by /u(Ar) dU ~ 1. 

The set of random unitary N x N matrices U with probability distribution given by the Haar measure dU, is called 
Dyson's Circular Unitary Ensemble (CUE) in N dimensions. The CUE plays a central role in the random-matrix 
modeling of chaotic and/or disordered systems with broken time-reversal symmetry. A new analytical method for 
treating the CUE is introduced below. For simplicity, the method will be illustrated at the example of the generating 
function 

F^,,r._m^})-l^^^dU n j3^^(^_^,,,„^)nDet(l-e---.t/t)' 

from which all information about the spectral correlations of the CUE can be extracted, if+a, '^-p and O-p 
are angular variables. (The latter are given a small imaginary part to prevent the determinants in the denominator 
from becoming singular.) The method does not exploit the invariance of i^„^,n_ under conjugation U i— » gUg~^, nor 
does it make any reference to the eigenbasis of U. For this reason it can be easily extended to the calculation of more 
general types of correlation function, such as those relating to wave functions and transport coefhcients. 

A. Generalized Hubbard-Stratonovitch Transformation 

The method we will use is a variant of Efetov's supersymmetry technique Q, which has become a standard tool 
for treating Hamiltonian systems with Gaussian disorder. Efetov's method starts by the usual trick of expressing the 
ratio of determinants as a Gaussian superintegral: 

F„^,„_({0,^}) 

dU j D{^,^;x,x) exp(-0tJ(5'^-'-e*''+°!7'=')</>V„-X+J<5^'-e'^+°t/'=')xU 

-<^^^(5" - e-'''-^C7"=)0L^ - X-/3('5''' - e-*'^-''C/"=)x'_;3), 



U(JV) 



where 1/)+^, 0+q,, and {k,l = 1,...,N; a — /3 — l,...,n_) are complex commuting (or bosonic) 

variables, with the bar denoting complex conjugation, and Xi X are anticommuting (or fermionic) variables. U'^'' are 
the matrix elements of the unitary N x N matrix U. The symbol D{(j), cj}; x, x) denotes the "flat" Berezin integration 
measure, i.e. it is given by the product of the differentials of the commuting variables times the product of the partial 
derivatives with respect to the anticommuting variables. The summation convention is used throughout this paper. 
The integral converges for 

Im^+Q > > Im6'_^. 

To simplify the notation we introduce a composite index a = (a, a) where a takes values = B (Bosons) or = F 
(Fermions). We then define a tensor ip with (n_ + n^)N —: nN components ip±a by 

ip is defined similarly. Setting D{(/), ip; x, x) = ^("^i '/') and combining the phases {9, ip} into a single symbol {lu} by 
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we may express the generating function in the abbreviated form 

Fn+,n. {{lo}) = [ Di4>, 4') I dU exp (-^ta(^'' - e^"+"C/")^Va 

J Jv(N) ^ 

-i'-biS''^ - e-'"-''C7"=)V'!_b) . (2) 

In the case of Hamiltonian systems, the next steps of Efetov's method are (i) to average over the Gaussian disorder, 
thereby producing a term quartic in the ip fields, (ii) to decouple this term by the introduction of a Hubbard- 
Stratonovich supermatrix field Q, and (iii) to make a saddle-point approximation eliminating the so-called massive 
modes. The final outcome of this convoluted and mathematically very delicate procedure is an integral over fields 
that live on a nonlinear supermanifold, called Efetov's cr-model space with unitary symmetry. It is difficult to see how 
this procedure could be adapted to the present case as the Haar integral over U(A^) is very far from being Gaussian. 
Fortunately, there is another way to proceed. 

The key innovation made here is to trade the integral over V{N) directly for an integral over Efetov's cr-model space 
with unitary symmetry. This is achieved by the following remarkable identity: 



V(N) 



dU exp (4>\JJ''4'la + ^i.t^^^V-- 
D^JiN{Z, Z) exp + i^UZba^ia) , (3) 

which can be proved by letting the Lie supergroup G\{nN\nN) act on a Fock space of bosons and fermions and using 
some basic concepts of the theory of generalized coherent states. The proof is given in detail in the Appendix. Here 
we confine ourselves to spelling out the meaning of the right-hand side. The new integration variables Z = {Zab} and 
Z = {Zab} are supermatrices whose Boson- Fermion block decomposition is written 



f Zbb 


Zbf \ 


z — ( 


Zbf \ 




Zff J 


\ Zfb 


Zff J 



Each block of Z (Z) is a rectangular matrix of dimension n+ x ri_ (n_ x The integration measure is 

Dfij^{Z, Z) = D{Z, Z) SDet(l - ZZ)^, 

where D{Z, Z) denotes the flat Berezin measure on the space of supermatrices Z, and SDet is the superdeterminant. 
The domain of integration is defined by the conditions 

Zbb ~ ^BBi Z-pF = ^^ffj 

and the requirement that all eigenvalues of the positive hermitian ri_ x n_ matrix ZbbZbb be less than unity. The 
integration measure is normalized by J Dfi^iZ, Z) = 1. In the next subsection we will argue that Z and Z should be 
interpreted as parameterizing Efetov's unitary cr-model space. 

Let me mention in passing that integrals of a type similar to the left-hand side of (^) appear in lattice gauge theory, 
where ip and U represent the quark and gluon fields, respectively. Motivated by this similarity, one might refer to 
the upper index of as "color" and the lower one as "flavor" . We could then say that (||) integrates out the gluon 
fields, which carry color, and replaces them by an integral over color-singlet meson fields carrying flavor. Note that 
the color degrees of freedom of the tensor ip are uncoupled on the right-hand side. 

By using the generalized Hubbard-Stratonovitch transformation (H) with i/);^^e'"+'' and "0^^ — > '!/'^i(,e~*"-'', 

we can process the expression (||) as follows: 



'^To be precise, I should say that this is true only locally. Berezin measures typically suffer from a global anomaly [ |ii| , [L2[ . 
What determines D{Z, Z) and its anomaly uniquely ( up to multiplication by a constant) is the requirement of invariance under 
translations on Efetov's cr-model space G/H, see Sec. [IB and the Appendix. 
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= j DfiN{Z,Z) SDct-^ 
= j D^iNiZ, Z) SDet"^ (l - Ze^^+Ze" 

D{Z, Z) SDet"^ f 1 - (1 - ZZ)-^Z{e"^+Ze~'"^- - Z)) . (4) 



In the first step we did the Gaussian superintegral over ip^ ip, producing a superdeterniinant. As the different species 
labeled by i — l,...,N are uncoupled, this superdeterniinant separates into a product of N factors. The expression 
e^'^+ Z stands for the supermatrix with matrix elements e^'^+'Zab- The second step in (Q) is immediate from the 
elementary identity 

SDet (^^ = SDet(A) SDet(L> - CA-^B). 

The third step follows from the formula D^fq{Z, Z) = D{Z, Z) SDet^(l — ZZ) along with the multiplicative property 
of the superdeterminant: SDet(A) SDet(B) = SDet(AB). 

Let me emphasize that the result (Q) is exact for all N . In contrast with Efetov's method, no saddle-point approxi- 
mation was used, and no massive modes had to be eliminated. Equation (Q) reduces the integral (|l]) over the iV^ real 
freedoms of U(A'^) to an integral over the 4n+n_ complex freedoms of the supermatrices Z, Z. Clearly, such a scheme 
is efficient when N is large and n+n_ is small. 



B. Geometric meaning of Z,Z 

We start with a concise mathematical description of Efetov's cr- model space of unitary symmetry. Let G = Gl(n|n), 
the complex Lie supergroup of regular supermatrices of dimension {n -\- n) x (n + n), and put n = ua + nji. (Later 
we shall make the identifications ua = and njj = n_.) Using tensor-product notation, we take G Gl(n|n) to 
be the diagonal matrix = ® diag(l,i^, — Imh) where is the n^i-dimensional unit matrix, and lx|i is the 
(1 + l)-dimensonal unit matrix in Boson-Fermion space. Matrices h G Gl(n|n) that commute with Sj. are of the form 

~^nR\nR J 

is the matrix presentation of S^. Thus the condition hY,z = Sz/i (or, equivalently, h = E^/iS^) fixes a subgroup 
H ~ G\{nA\nA) x Gl(nfi\nji) of G. The coset space G/H is a complex-analytic supermanifold in the sense of Berezin- 
Kostant-Leites [ p^ , p^ , and has complex dimension dBoson = c^Fermion = '^nAriR. The canonical projection G G/H 
endows the coset space G/H with a natural G-invariant geometry. If G/H is modeled by supermatrices Q = g'Szg"^, 
this geometry is given by the rank-two tensor field STr(d(5)^, where STr denotes the supertrace. The G-invariant 
Berezin integration measure on G/H is denoted by Dgn, or DQ. The objects one wants to integrate are the 
holomorphic functions on G/H, i.e. functions with a holomorphic dependence on a set of complex local coordinates 
of the complex-analytic supermanifold G/H. Such functions are written f{gH) or F{Q). The domain of integration 
is taken to be a Riemannian submanifold Mb x M-p of the support of G/H , where 

Mb = V{nA,nR)/V{nA) x U(nflJ (BB sector), 
Mf = lJ{nA + nflJ/U(nA) x U(nfl,) (FF sector). 

In this way one gets a G-invariant Berezin integral 

/ Dgn f{gH) = f Dgn figogH) {go G G), (5) 

J Mb X Mf J Mb x Mf 

which is called the integral over Efetov's a-model space with unitary symmetry. Translation gH l—^ gogH by an element 
go € G does not leave the integration domain Mb x Mp invariant, in general. Nevertheless, if / is holomorphic, as 
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is assumed, (||) does hold because Cauchy's theorem appHes and allows to undo the deformation of the integration 
domain. In the terminology of the pair {G/H,AIb x Mp) is called a Riemannian symmetric superspace of type 
^III|^III. (Nine more types e xist. ) 



To make contact with Sec. [I A we introduce a suitable parameterization of G/H . Using the above presentation 



where = diag(l„^|„^, -1„„|„^) we decompose g e G as 

f A B 

^=\C D 

Since h — diag(/i+, the invariants for the right action of ft. £ iJ on G are Z := BD~^ and Z :— CA~^. We may 
take the matrix elements of Z, Z for a set of complex local coordinates of the coset space G/H. The expression for 
Q — gYizg~^ in these coordinates is 

1 Z\ f 1 \ / 1 Z^ 



^={z -i){z i) ■ 

From the expression for the invariant rank-two tensor field, 

STr(dQ)2 = const x STr(l - ZZ)-^dZ{l - ZZy^dZ, 

one easily finds the G-invariant Berezin measure DQ in these coordinates to be the (locally) flat one, D{Z,Z) 
The flatness results from cancellations due to supersymmetry. 

To make the restriction to the Riemannian submanifold Mb x M-p, note that for a unitary matrix 9 — d) ^ 
V{nA + rifi), we have 6^ — — ((f — ca^^b)^^ca^^ and d'^ — {d ~ ca^^b)^^. This implies 

zip = (M"^)t = d-^h^ = -ca-^ = -Zff- 
For the noncompact case g G \J{nA, tir), the expression for fe^ has its sign reversed, so 



^BB = ^"^^^^ = = +Z: 



BB- 



Moreover, the pseudo-unitarity of g G \J{nA,nii) implies that a^a — c^c — 1, from which we infer 

ZbbZbb = ca^^ a^^'^ = c(l + c^c)^^c\ 

so that all eigenvalues of the positive hermitian ur x ur matrix ZbbZbb must be less than unity. 

All these facts suggest to put ua = tir = n_ and identify the supermatrices Z,Z defined here with those 
figuring in the generalized Hubbard-Stratonovitch transformation (^). This identification is made rigorous in the 
Appendix, where a Lie-algebraic proof of (^ is given. The proof relies in an essential way on the existence of the 
invariant integral (^, and on the identification of D{Z, Z) with the G-invariant Berezin measure Dgjj on G/H. 

For further orientation, consider the simplest example n+ = — 1, where Mb — U(l,l)/U(l) x U(l) ~ 
(two-hyperboloid) and Mp = U(2)/U(l) x U(l) ~ S^ (two-sphere). In this case the BB and FF blocks of Z are 
just numbers. The expression for Zff in terms of the usual polar (9) and azimuthal (ip) angles on S^ turns out 
[^6| to be ZpF — tan{9/2) expiip. From this it is easily seen that Zff can be interpreted as being a complex 
stereographic coordinate for S^ . Its BB analog can be parameterized in terms of the hyperbolic polar angle 9b on H'^ 

by ^BB = tanh(6'B/2) expi(/3B- 

To conclude this subsection, let us recast the result (||) in the coordinate-free language developed here. From (|^) or 

Z\f {1-ZZ)-' -Z{l-ZZ)-'\ _JQ++ 
^-3^^9 -[Z l)\z{l-ZZ)-^ -{l-ZZ)-^ ) ^■\Q-+ 

we read off the relations 



z{i - zzr^ = Q_+/2, (1 - zzr^ = (1 + g++)/2, 
z = -{i + Q++)-ig+_, (1 - zz)-^ = (1 - Q—)/2. 

Inserting these into (jj) we can express the generating function as a Q-integral: 
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^;+,„_(M)= y Z?OSDct-^(l-0__+Q_+e''^+(l + Q++)-ig+_e-''^-). (7) 

Let me repeat that this expression is exact for all N. 

The following remark may be helpful. It is well known jl^ that the density of states of the Gaussian Unitary 
Ensemble tends to a semicircle in the large- limit. It is also known jl8| that the so-called "massive modes", 
which show up in Efetov's method and are eliminated by a saddle-point approximation, are needed to reproduce the 
nonuniformity of the semicircular density of states. In contradistinction, the eigenvalues of a unitary matrix U S CUE 
are uniformly distributed over the unit circle, on average. Qualitatively speaking, it is this uniformity of the CUE 
spectrum that allows a Q-integral representation to be derived without introducing the massive modes and without 
making any saddle-point approximation. 



C. Large-N Limit 

For further calculations the formula (^) turns out to be most convenient. In the large- limit the integrand in (Q) 
substantially differs from zero only when the difference of matrices uj+ — cj_ is of order This fact allows us to 

expand: 

e"^+Ze-"^- -Z =^i {uj+Z - Zuj^) + ©(l/A^^), 
which leads to the approximate equality 

Fn^^nAi'A)- Jd{Z,Z) expiNSTr(uj+ZZ{l- ZZy^ -u^ZZ{1- ZZy 

The N eigenvalues of the unitary matrix U are distributed over the unit circle with circumference 27r, which results 
in the mean density of eigenphases being ly = N/2-k. By setting uj :— diag(a;+,w_) we can write the integrand in the 
form expiiVSTra;((3 — Sz)/2. Thus, the large-A^ limit of the generating function F„^_„_(ijj = {oj}) evaluated on a 
scale set by the mean spacing, is 

lim F„ „_ (27rti;/A'^) = / dQ exp z7rSTrcL'((5 — S^). 

Equivalently, we may write Q ~ gYi^g^^, and get 

lim F„^,„_(2^c^/7V) = / i?gffexpi^STrtj(5E,g-i -S,). (8) 

As is well known [p| Jl8|jT^ , the same expression is obtained for the Gaussian Unitary Ensemble, which confirms the 
expected equivalence |10| of the Circular and Gaussian ensembles in the large- limit. 

We are now going to calculate from the large- result (^) the so-called rt-level correlation function i?„ for any n. 
(Efetov H] calculated only the two- level correlation function.) To that end, we will use a supersymmetric variant of 
a result known in the physics literature as the Itzykson-Zuber formula ]l9| , which is actually a special case of a more 
general result proved 23 years earlier by the mathematician Harish-Chandra. 

Let K be any connected compact semisimple Lie group, which we assume to be realized by n x n matrices k, and let 
T be a maximal commuting (or Cartan) subalgebra of Lie(i^). Without loss we take T to be the diagonal matrices 
in lAe{K). If A, B are two elements of T, Harish-Chandra showed that (see Theorem 2 of pfl]) 



K 



dkexpTvAkBk^ = ——— V (-l)!"! exp TrAsB, (9) 
p(A)p(B) ^ ' ' ^ ^ 



where p{A) — Y\a>o'^i^) product over a set of positive roots a of the pair [T, Lie(i4r)], and const is a 

normalization constant. The sum runs over the Weyl group W[ivr], which is the discrete group of inequivalent 
transformations T — > T by A i-^ sA := sAs~^ with s € K. The Weyl group can be considered as a subgroup of the 
symmetric group S„ acting on the elements of the diagonal matrices A d T. The symbol \s\ = 0, 1 denotes the parity 
of s. For the special case K = SU(n) (or U(n), it makes no real difference) one has p{A) = Y[i<ji^i ~ where 
Ai are the elements of the diagonal matrix A, and the Weyl group coincides with S„, i.e. {sA)i = {sAs~^)i =: ^^(i) 
where s G S„ is a permutation and |s| is the parity of that permutation. The general result (^ thus takes the more 
familiar form ||l^ 
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/ dUcxpTrAUBU^^ = const x T\{Ai - Aj)-^{B^ - Bj)-^ x Det (e^'^^j . .^^ 



The right-hand side of (|) can be viewed as the stationary-phase approximation to the integral on the left-hand side, 
with the points of stationarity being enumerated by the elements s G W[K]. Thus, Harish-Chandra's formula states 
that the stationary-phase approximation is exact in this case. M. Stone has traced this remarkable coincidence 
to a hidden symmetry of the integral. 

I conjecture that the formula (S) extends to any compact classical Lie supergroup K: 

Dk expSTiAkBk-^ = ——— ^ (-1)1^1 exp STrAsB, (10) 

■^^ ' seWlKo] 

where W[ifo] now is the Weyl group of the ordinary group Kq supporting K, and p{A) turns into the product of 
positive bosonic roots divided by the product of positive fcrmionic roots. The strategy of the proof ought to be a 
supersymmetric extension of that of Harish-Chandra and use the theory of invariant differential operators ||20| , ^ , ^ . 
For the special case of the unitary Lie supergroup K — U(n|n), which is the one we will need below, two proofs using 
alternative strategies can be found in Refs. and |25|. In that case, if A = diag(>li_B, A„^b; ^i,f, ■•■,^n,F), one 
has 

p{A) = n(A,,B - A,-b)(A,,f - A,,p)/Y[{A,B - A,-f), (11) 

and W[Ko] = W[U(n) x U(n)] ~ S„ x S„, so 

V (-1)1^1 expSTrAsB = Det (e^''^-^^-^) . . , x Det (e"^' '''^^ '') . . , 

Z ^ \ / x- \ / t,j — l,...,n \ — l,...,n 

seW[Ko] 

In a reduced form tailored to physical applications, the Harish-Chandra-Itzykson-Zuber integral for U(n|n) made its 
first appearance in |2^ . 

Now recall that the n-level correlation function i?„(0i,6'2, ...,6'„) is defined as the probability density to find, 
given a level (i.e. an eigenphase) at 0i, n — 1 levels at positions 02, ■■■,0m irrespective of the positions of all other 
levels. By making use of the identities 



d_ 

dip 



H4l:i^..Tr(l-e-C/t)-\ 
¥.=eDet(l -e*»C/) ^ ' 

~ 27-77:7 — ' '^(^)' 



(1 - e-^-»»)(l - e-^+»») 2(1 - cos6l) -h 4sinh^(e/2) 
we can extract Rn from the generating function ivi,„, Eq. (|l|), by the following formula: 

n 

i?„(6'i, ...,6i„) = lim - Fn,n{{9a+ie,9a~ie,p+a,'P-a}), (12) 

which serves as the starting point of our supersymmetric calculation of i?„. 

The first step toward doing the integral (§|) is to introduce Efetov's polar coordinates |^ on the coset space G/H: 

gH = haH, where ^ ^ 'q^ h ^ ^ 

and a = exp^^ z = diag(a;i, a;„, zyi, iy„). 

The variables Xi,yj contained in a are the "radial" coordinates, and the supermatrix h is parameterized by a set of 
"angular" coordinates. The choice of integration domain Mb x Mp for G/H dictates that we must integrate over real 
Xi,yj and superunitary h G U(n|n) x U(n|n) C H. The polar-coordinate form of the invariant integral (||) is 



jDgHfigH) = j (^j Dh f{haH)^ J{a)da + ..., 



(13) 



where Dh is the Haar-Berezin measure of H, and da denotes the Euchdean measure on the abehan group of radial 
elements a. The Jacobian of the transformation to polar coordinates is given (2^j2^ by J{a)da = 



U^^j smh {xi-Xj) smh {x, + Xj)sm {y,-yj)siii (yi + y^) ^ 

— —27 , . X . , 2/ — ^ 11 smh(2a;fc) siii[2yk)dxkdyk 

Ihj smh [x^ + zy^) smh (xi - lyj) 

Y\^^Acofi\\2x^ — cosh2a;j)^(cos2j/i — cos2yj)'^ 

= const X i = — — d(cosh2xfc)d(cos2?/fc). (14) 

J (cosh 2xi - cos 2yjy j^J-^ 

For the second equality sign the trigonometric identity 

2 sinh(Q! + (3) sinh(Q! — /3) = cosh 2a — cosh 2/3 

was used. The integration domain is taken to be < < oo, < j/fc < 7r/2. The dots in ( p^ indicate correction terms 
that arise from the anomaly of the invariant Berezin measure in polar coordinates. The form of these corrections was 
investigated and completely determined in [ p9| . Their main characteristic is that they are supported on the boundary 
of the integration domain (which is why they are alternatively known as "boundary terms"). More precisely, they are 
obtained by setting Xi = yj = for one or several pairs {xi,yj). It will turn out that these boundary terms do not 
contribute to Eq. ([l^), and therefore do not need to be specified explicitly for our purposes. 

The next step is to express STrcjgE^g^^ in the exponential of the integrand of (^ in terms of polar coordinates: 

STTugY^zg^^ = STrT, z^ha'^h^^ =^ STrE^cj/i cosh(2 In a)ft,^"'^ 

= STrLu^h^Ah:^^ - STrLu_h_AhZ^, 

where A = cosh22;, lu = diag(cij+, lu-), and uj± = diag(0i ± ie, 0„ ± is; </?±i, </3±n)- By this and Dh ~ Dh^Dh^ 
the angular integral over U(n|n) x U(n|n) factors into a product of two integrals over U(n|n). By formula (^0|) for 
\]{n\n), the integral over gives 

/ D/i+expiSTrw, /i, ^/i-i = — _— y"(-l)l^"l exp iSTrtj+M. 

J\J{n\n) P{^+)P{A) ^ 

A similar expression results from doing the integral over Multiplying both integrals we get a factor p{A)~'^ . 
By (|ll) and A — diag(cosh2a::i, .. cosh2x„;cos2j7i, ...,cos2jy„) this factor exactly cancels the factor multiplying 
dA :— d(cosh 2xfc)d(cos 2yfc) in (|lj). Hence we obtain 

DgH cxp iSTrtjgEzg"^ 
const 



p{uj+)p{uj-) ^ 



boundary terms. 



The inverse of the product of positive roots p{uj+)^^ vanishes linearly with each difference w+a,B ~ ^jJ+a.F = Oa + is — 
tp+g. An analogous statement holds for p(cj_)^^. Therefore, on taking the 2n first derivatives at (p±a — 9a±is [see 
(|l2|)], we simply obtain a constant, which can be absorbed into the normalization. Note 'STT:YjzLo\^^^^Q^±i^ — 0. 

A further simplification results from the fact that dA and the integration domain for A are Weyl-invariant. This 
permits us to reduce the double sum over s+, s_ to a single sum over the relative permutation s sZ^s^, times an 
overall factor dim W[U(n) x U(n)] = dim(S„ x S„) = rt!^. Absorbing this constant factor into the normalization, we 
obtain for the n-level correlation function the following intermediate result: 



Rn = const X ^(-1)1^1 Vans'" j dAexpiSTrA {{6 + is) - s{ 



IS, 



where = diag(6'B, ^f), O-q = 0-p = diag(0i, 6'„). The integral factors into a product of 2n one-dimensional 
exponential integrals over the elements of the diagonal matrix A, n of which are compact (—1 < cos2yi < 1), and the 
other n are noncompact (1 < cosh2a::i < oo). The Weyl group element s =: (sb,Sf) £ W[U(n) x U(n)] acts on the 
BB sector by sb and on the FF sector by sp. If 9i = (sB0)i, the integral over ti = cosh2a;i diverges as s~^ in the limit 
s ^ 0. The maximal divergence occurs for sb = identity, in which case a singular factor e^" is produced, canceling 
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the prefactor e" and producing a finite result in the limit e 0. The boundary terms mentioned earlier disappear at 
this stage as they contain at most n — 1 noncompact integrations and therefore diverge more weakly than What 
remains are the n compact integrations over ti = cos 2yi. Using the elementary integral J^^ dt exp iirOt = 2 sin(7r0) / 7:9 
we immediately arrive at the final result: 

seW[U(n)] 1=1 ^ ^ ^> \ \ I JJ y i,j = l....,n 

where the correct value of the normalization constant was restored by hand. This result coincides with the one 



obtained by the Dyson-Mehta orthogonal polynomial method ||17[ in the large- TV limit. In Sec. |V| we are going to see 
that the first form of the result ( p^ ) extends in a very simple way to the circular ensemble that is obtained by taking 
instead of the unitary group the symplectic one. In that sense, this is the "good" way of writing the result for 



III. CIRCULAR ORTHOGONAL ENSEMBLE 



If an open quantum mechanical system possesses an anti-unitary symmetry T (time reversal, for example) and 
= +1, there exists a basis of scattering states such that the scattering matrix 5* is symmetric: S = S'^ . The set of 
symmetric S'- matrices can be parameterized in terms of the unitary matrices U £ V{N) by S' = UU^ . The product 
S = UU'^ is invariant under right multiplication of U by any orthogonal matrix, which means that S = UU'^ 
lives on the coset space U(A^)/0(A^). 

The Circular Orthogonal Ensemble (COE) of random-matrix theory is defined by taking S to be distributed 
according to the uniform measure on U(iV)/0(A^). We denote this measure by dfi{S). The generating function for 
COE spectral correlators, Fn^^n^, is defined as in (0) but with the substitutions U(iV) -> U(7V)/0(iV), U ^ S, 
and dU — > d^{S). Because the invariant measure dfi{S) is induced by the Haar measure dU of U{N) through the 
projection U{A^) \J{N)/0{N), we can write the defining expression for J^„_^,„_ also as follows: 

F ((f) <n^^ [ dJTT\ Det(l-e-'^^°t/[/^) Bcijl - e~^^->'^UW) 
{{e, ^}) - j^^^^ dU [[ ^-^^—^-^ Det(l-e-^-.f7m)- 

This expression will now be processed by the method of Sec. ||. To do so we use the trick of doubling the dimension: 

Det (1 - e^UU'^) = Dct 

The extra degree of freedom implied by this doubling will be called "quasispin" and denoted by the symbols t, \.. 

After the introduction of quasispin, we express i^n_|_,„_ as a Gaussian superintegral in the usual way, see Sec. [I A 
Eq. (^. As before we put {w} = {9, ip}. Then we lump the two terms containing U into a single one: 

where the tensors 0+ and x+ have components 

Here |a| = if a = (a, B), and |a| = 1 if a = (a, F). The terms involving U'^ are manipulated in a similar manner: 

+ ^L,^e-^^-'{U^y^^U^ ^ X^bU^'^-b, where 

In the next step we apply the generalized Hubbard-Stratonovitch transformation (^: 

J dUexp(^4>\AU''x\A + X^BU''cf,^-B) = J DtiN{Z,Z)exp(^cf,^_^^Z^s4>-B+X^BZBAxU 

Note that the exponential on the right-hand side separates into terms containing and x- The same is true for 
the terms that do not involve U. Therefore the integral splits into two integrals, one over the tensor ip and the 
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other one over the tensor x- Consider the integral first, and define an orthogonal matrix t+ ("charge conjugation") 
by requiring the tensor (t+)^^,(/>!,_^, to have components {^!|-a-f ) V'+oT }■ The matrix that does this job is given in 
quasispin block decomposition by 



T+ = 



a 

1 



where a is the matrix for superparity {c7aa' = (— l)'"'^aa')- The same matrix acting in the space of "negative" indices 
(with the appropriate change of dimension) is denoted by r_ . Having introduced these matrices we can write 

i}\a^^\a^ = |'^Va(^+)aa'<^+A'' and V;!_6|V-6T = 5'/'-b(^-)bs'<^-b'- 
On collecting terms we encounter the integral 

j Dcf) exp {4^+aZab4'-B - W+A{r+)AA'4^+A' - W-B{r-)BB'4'-B') 

= (-l)"^/2SDet-^/2 ( _ J+T^ = SDet-^/^ (l - Z^r+ZTr') • 

Here {Z'^)ba = ^ab(-1)^'^'+^^'^' is the supertranspose, with |^| = if A = +(a,B) t, and |A| = 1 if A = +(a,F) t- 
The X integral is done in a similar fashion and gives SDet-^/2 (l - ZT^^e^'^+Z'^e-^'^-T_y By making the variable 
substitution 

Z ^ e*"+/^iZe-*'^-/'^, Z ^ -e''^-^SZe''^+/'^, 
we bring the final result into the symmetrical form 

F„^,nAW})= j Dm{Z,Z) SDet-^/2 (i + z'^T+e'''+/^Ze-"^-/^Tl^) 

xSDet-^/2 ^ ^^-lgia.+ /2^Tg-»a;_/2^_j ^ 

which is exact for all A''. The matrices Z,Z parameterize the complex coset space Gl(2n|2n)/ Gl(2n+|2n+) x 

Gl(2n_|2n_). 

We turn to the large- iV limit. In order for this limit to be nontrivial, we must again take the difi'erence lu^ — ui- to 
be of order 1/N. By Dfi^iZ, Z) = D{Z, Z) exp A''STrln(l — ZZ) the integrand at a;+ = a;_ = const x identity can be 
written as 

exp (^AfSTrln(l - ZZ) - ySTrln(l + Z'^t+ZtZ^) - ySTrln(l + Zt+^Z'^t_)^ . 
For large A'', the integral is dominated by contributions from the subspace determined by the equation 



1 vT, 



Z = -T^'Z 



This subspace is referred to as the saddle-point manifold. Integration over the Gaussian fiuctuations normal to the 
saddle-point manifold gives unity by supersymmetry, in the large- iV limit. 

By restricting the integrand to the saddle-point manifold, expanding with respect to uj+ — co- and keeping only the 
terms that survive for A'' ^ oo, we obtain 

^lim F„_^,n_(27r{u;}/iV) = j D{Z, Z) expinSTr (uj+ZZ{l - ZZ)-^ - lo_ZZ{1 - ZZy'^y 

This has the same form as for the CUE, except for the replacement N N/2, the doubling of the dimensions of the 
supermatrices Z,Z, and the imposition of the constraint Z = —t^^Z'^t_. On putting cj = {u>} = diag(t<;+, W-) we 
can re-express the result in the form 



lim F„+,„_(27ra;/7V) = / Dgn expi7rSTiu{gi:,g-^ - 



)/2. 
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The relation between Q = gT,zg ^ and Z,Z is formally the same as before. The constraint Z = —t_^_^ Z'^t_ is 
known ||l^,|l^ to reduce the complex superspace Gl(2n|2n)/Gl(2n+|2n+) x Gl(2n_|2n_) to the complex submanifold 
Osp(2n|2n)/Osp(2n+|2n+) x Osp(2n_ |2n_). The conjugation relations Zgg = Zgg and Zpp = —Z^p translate into 

Mb = SO(2n+,2n_)/SO(2n+) x S0(2n_), Afp = Sp(2n+ + 2n_)/Sp(2n+) x Sp(2n_). 

For — ri- = 1, the above expression for the generating function Fn_^^n- can be calculated by introducing Efetov's 
polar coordinates. The resulting formula for the two- level correlation function R2 coincides with that of ||^. It is not 
clear at present how to do this calculation for general n+, rt_. [Harish-Chandra's formula does not extend to this case 
as neither nor lo = diag(ti;+, w-) are elements of Lie (Osp(2ri|2n)).] 



IV. CIRCULAR ENSEMBLE OF TYPE C 



Dyson's Circular Ensembles COE and CSE are constructed by starting from the unitary group and passing to the 
coset spaces \J{N)/0{N) and U(2A^)/Sp(2A^). Other circular ensembles of relevance ||31] to mcsoscopic physics can 
be obtained by taking for the starting point the orthogonal or symplectic group instead of the unitary one. Again 
the option of forming various coset spaces exists. Of these possibilities the technically simplest one is to consider 
an ensemble of scattering matrices S drawn at random from the symplectic group Sp(2A^), with no coset projection 
done. This is what we shall do in the present section. The probability distribution for S will be taken to be the 
Haar measure dS of Sp(2iV). The circular ensemble so defined is called "type C". A physical system where such 
an ensemble can be realized is a chaotic Andreev quantum dot ||3^ with time-reversal symmetry broken by a weak 
magnetic field. 

The defining equations of Sp(2A^) are 

S-^^ = S = CS-^'^C-\ where C = iay(E)lN=( ? 

is the symplectic unit in 2A^ dimensions. The polar (or Cartan) decomposition of an element S £ Sp(2A^) has the 

form S — ke^^k~^ where 9 — az® diag((?i, 6'jv) and cr^ — diag(-|-l, —1). Thus if 9i is an eigenphase of 5, then so is 
—9i. This symmetry under reflection 9 —9 is called a "particle-hole" symmetry. 

All information about the eigenphase correlations of the circular random-matrix ensemble of type C is contained 
in the generating function 

isp(2JV) ^J-^ Dct(l -e^'^-S') 

This expression is simpler than the corresponding one for the CUE because Dct(l — zS") = Det(l — zS't) as a result of the 
particle-hole symmetry. By putting {uia\ = {Scn'^a} and proceeding as in Sec. II A| , we get a Gaussian superintegral 
representation for Fn'. 



F^{{uj]) = / / dS exp-^^(5'J - e'^-U'^Wa. 

J JSp(2N} 

The components ipl^, "ipl^ of the tensors ip, are bosonic for a = (a, B) and fermionic for a — (a, F). The analog of (|^) 
reads: 

d^exp^^^'^^^ - / i?MA^(Z, Z) exp [W^Z^^C^ - Wai^lt^ Z ^^C^ . 

Sp(2Af) JMbxAIf ^ ' 

As before, DfiN{Z, Z) = D{Z, Z) SDet(l — ZZ)^ , but Z and Z are now square supermatrices and are subject to the 
symmetry conditions 

Zab = (-l)l''l + l^l + l'^ll^l+lZb,, Zab = (-l)l'^ll''l + lZb,. 

These conditions can be succinctly written as 

Z = -Z'^a, Z = -aZ^ , 
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where a is the superparity as before. They express the fact that Z, Z parameterize the complex supermanifold 
G/H — Osp(2n|2n)/Gl(n|n), see also below. The integration domain Mb x Mp is fixed by 

^BB ^ ^BB' ^FF ^ ^FF' 

and the requirement that all eigenvalues of the positive hermitian n x n matrix ZbbZbb be less than unity. It can be 
shown |l^ that this means 

Mb = S0*(2n)/U(n), Mp = Sp(2n)/U(n). 

{G/H, Mb X Mp) is called a Riemannian symmetric superspace of type DIIIICI. 
In the simplest case n = I the supermatrices Z, Z have the form 

z=[? ^A, z=( 



where 21,^2 are complex commuting and CijC2 are complex anticommuting numbers. It is seen that the BB sector 
Mb is trivial in this case (in fact S0*(2)/U(l) consists of just a single point). The variable Z2 = ~z\ of the FF sector 
can be interpreted as the complex stereographic coordinate of a two-sphere ~ Sp(2)/U(l). It is very important to 
note |]l2| that the variables ^1 , C2 are not related by any kind of complex conjugation. 

The proof of the generalized Hubbard-Stratonovitch identity for the present case is closely analogous to that for 
the CUE, presented in detail in the Appendix. (To understand the following remark you must study that appendix 
first.) The key idea is to consider the generalized coherent states 

\N/2 



\Z) := exp (\-e,Z,,C^-c^,) |0) SDet(l - ZZy 
{Z\ := SDct(l - ZZ)^/2 (0| exp (i<(-l)l'^IZ,,C^V,) , 



built from an absolute vacuum cjj|0) = by repeatedly acting with Bose and Fermi creation and annihilation operators 
obeying the supercommutation relations [c^,c^] = S'^^Sab, and to exploit the fact that P = J D{Z, Z)\Z){Z\ projects 
on the singlet sector of the symplectic group acting on Fock space by 1-^ S^^c{. 
Use of the generalized Hubbard-Stratonovitch transformation leads to 



1 -e^^Ze'- 
-Z 1 



Fn{M) = J D^iN{Z,Z) SDct- 

= j D{Z,Z) SDet-^ " (1 " ZZy^Zie'^^Ze"^ - Z)j , 

where uj stands for the diagonal matrix with elements uja- To arrive at an equivalent Q-integral representation, we let 
g run through the orthosymplectic Lie supergroup G — Osp(2n|2n) defined by 



^-^GD -[aO [B^D^I I 1 



and set Z = BD-^, Z = GA'^ 



I Z\ 1 W 1 Z\ Q++ g^ 



z I ) \ Q -I J \z I J • VQ-+ Q 



as before. The subgroup H d G oi elements h — diag(A,v4^^ ) satisfying hY^zh^^ = '^z is isomorphic to G\{n\n). 
Hence the elements of the supermatrix Q are functions on the coset space G/H = Osp(2r7,|2n)/Gl(n|n). By simple 
manipulations we obtain 

K(M) = /^Q SDet-^ (1 - Q__ + Q-+e'^{l + Q++)-'Q+-e'") , 

which is to be compared with (|^). 

In the large- limit the exact expression for F„ simplifies to 
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lim Fn{TT{uj}/N) = [ DgHexpiTTSTTu;igj:,g-^ - E,)/2, 

where w = diag(w, —uj). Note that this is formally identical to (ph. The CUE n-level correlation function i?„ was 



extracted from that equation in Sec. [IC. It turns out that the corresponding calculation can be done in the present 
case, too, by extending Harish-Chandra's formula to the orthosymplectic Lie supergroup. [In contrast with the 
situation for the COE, the matrices lo and now are elements of Lie (Osp(2n|2n)).] To keep th e pre sent paper within 
size, let me omit the details of this calculation (which in any case is very similar to that of Sec. |ll C| ) and simply state 
the result for R„ thus obtained: 



lim Rninei/N, 7r0„/A^) - V (-l)''' [T 

seW[Sp(2n)] i=l 



The Weyl group W[Sp(2n)] is generated by the operations of transposition 9i-i ^ = 2,...,rt) and reflection 

9i 1-^ —Oi, all of which have odd parity \s\ = 1. Please observe the perfect analogy to (Iq): the Weyl group of U(n) 
has simply been replaced by the Weyl group of Sp(2n). Specialization to n = 1, 2 gives 



sin 2'k9 



1 



hm R,{7:9,/N) = l- , 
N^oo Znui 

^ , r. i^r r. / , rx sin27r6'i\ /, sin27r6'2 
lim R2{tt9i/N,tt92/N) =1 — ^ 1 



'sin7r(6'i - 6*2) sin tt (6*1 + 6*2) 
n{9i - 92) n{9i+92) 

Precisely the same expressions were obtained for the correlation functions of the Gaussian random-matrix ensemble 
of type C in pl[ | . The reasoning used there was heuristic and took recourse to the mapping on a one-dimcnsional free 
Fermi gas with Dirichlet boundary conditions at the origin 9 = 0. This coincidence of results is of course expected 
from the large- equivalence of the Circular and Gaussian Ensembles. 



V. SUMMARY AND OUTLOOK 

The key result of this paper is the Hubbard-Stratonovitch identity (|^), relating an integral over the unitary group 
U(-/V) to an integral over a Riemannian symmetric superspace of type ^III|^III (Efetov's unitary cr-model space). The 
detailed proof given in the Appendix derives this identity from the standard properties of generalized coherent states. 
Instead of U(iV) one can also consider the symplectic group Sp(2A^), or the orthogonal group in an even number of 
dimensions, S0(2iV). In these cases, identities similar to (^) exist and relate the corresponding group integrals to 
integrals over Riemannian symmetric superspaces of type DllljCI and CI|-DIII respectively. For Sp(2A^) this was 



described briefly in Sec. IV. The case of S0(2A^) was not discussed here and is left for future work. 

In the present paper the generalized Hubbard-Stratonovitch transformation was applied to Dyson's Circular Unitary 
Ensemble (CUE, or type A) and the circular ensemble of type C. In both cases the large- limit of the n-level 
correlation function was calculated for all n, by appropriate extensions of Harish-Chandra's formula. The method 
is not restricted to spectral correlations but can be used for wave amplitude correlations and transport coefficients 
as well. It can also be adapted to other symmetry classes. The trick that works for Dyson's Circular Orthogonal 
Ensemble (COE, or type Al) is to write the elements S of the COE in the form S = UU^ with U £ CUE and 
then proceed as in the unitary case. A similar trick works for Dyson's Circular Symplectic Ensemble (type All), the 
Circular Ensemble of type CI and, in fact, any one of the large class ||lj] of circular ensembles. 

As was mentioned in the introduction, nontrivial applications to random network models are expected. What puts 
these within reach is the fact that the Hubbard-Stratonovitch scheme (||) is valid for all N, including = 1. This 
permits to transform a network model with random U(l) phases into a theory of coupled supermatrices Z, Z. Doing so 
for the Chalker-Coddington model ||^, for example, and taking a continuum limit, one obtains the two-dimensional 
nonlinear a model augmented by Pruisken's topological term. 

Another realm of fruitful application of the identity (^) will be quantum chaotic maps. Andreev et al. [Q have 
recently argued that energy averaging for deterministic Hamiltonian systems can be used to derive a "ballistic" |^ 
nonlinear a model, in much the same way as impurity averaging for disordered systems leads to the usual "diffusive" 
a model. When one goes from Hamiltonian systems to periodically driven ones or maps, the role of the Hamiltonian 
passes to a unitary operator U, the Floquet or time-evolution operator. As the spectrum of U lies on the unit circle in 
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C, an approach analogous to that of Andreev et al. will employ U(l)-phase averaging instead of energy averaging. It 
is clear that, by carrying out this phase average with the help of the identity (^), one will be able to derive a nonlinear 
(7 model for maps. Work in this direction is in progress |^^. 

Finally, let me mention that the identity (|^) is not restricted to the supcrsymmetric case. Similar formulas can 
also be derived rigorously when the spinor ip is purely bosonic or fcrmionic. (In the bosonic case, convergence of 
all integrals places a lower bound on the allowed values of N .) The fermionic version of is particularly exciting 
as it promises nontrivial applications to lattice gauge theory, by enabling an exact transformation from the gauge 
degrees of freedom to meson fields. (Note that the usual plaquette term TiUUUU of lattice gauge theory can be 
generated by coupling a number of heavy fermionic ghosts to the gauge field via ^pUijj.) In this way, one may succeed 
in constructing a nonperturbative proof of a result due to Witten |34[| , who summed planar diagrams to argue that 
large- iV quantum chromodynamics is equivalent to a weakly coupled theory of mesons, in the low-energy limit. Sadly, 
although the exact gluon-meson transformation works for the gauge groups U(7V) and Sp(2iV) (and, presumably, for 
S0(2iV) as well), it does not easily extend to SU(iV), except for N = 2, where one has the accidental isomorphism 
SU(2) ~ Sp(2). 

This work was supported in part by the National Science Foundation under Grant No. PHY94-07194, and by the 
Deutsche Forschungsgemeinschaft, SFB 341, Koln- Aachen- Jiilich. 



VI. APPENDIX 



Given a compact Lie group K and a unitary irreducible action T of K on a module V with lowest (or highest) 
weight |0), one can consider objects of the form Tk\0) {k G K). These are called generalized coherent states, and they 
have many nice properties Ka-B^. Most importantly, they resolve the unit operator on V: 



iv^ / dfcr,|o)(o|r2. 



K 

2 



where dk denotes the Haar measure of normalized by dfc|(0|Tfc|0)| = 1. A well-known application of this 
is the derivation of coherent-state path integrals for quantum spin systems, in which case K = SU(2) and y is a 
spin-S" representation space of SU(2). The objective of the present appendix is to extend some of the mathematics of 
generalized coherent states to the case where instead of K we have the supergroup G — Gl(n|n) that acts on Efetov's 
complexified coset space G/H of unitary symmetry. This will eventually yield an easy proof of the identity (|3|). 
We start by defining a group action of Gl(n|n) on a Fock space of bosons and fermions and develop ing variou s useful 



structures that come with it. The general context and our notations have been laid down in Sees. II A, II B, Recall 
that the tensor ?/; has components V'+aj V-'lfci where i = 1, A'^, and a = (a, cr), h = (/?, a) are composite indices with 
range a — l,...,7i_|_; [3 — l,...,n_; and cr = B,F (Bosons and Fermions). For notational convenience, we introduce 
a composite index A taking values +a or —6. Thus, {^Z"^} = {V'+a; ''/'-bl- The terminology we use is to refer to the 
upper index i as "color" and the lower index A as "flavor" . Now we introduce creation and annihilation operators 
c^,c^ for bosons and fermions, which obey the canonical supercommutation relations 

where |A| = if A = ±(q;,B) and |^| = 1 if A = ±(q;,F). These operators act on a Fock space with vacuum |0) 
defined by c^|0) = 0. They can be viewed as the quantum mechanical counterparts of the "classical" variables 'tl)\,^\. 
For the following it is very convenient (though not necessary) to make the transformation 

which is canonical, i.e. leaves the supercommutation relations unchanged. (The operators c!|_q, c!|„q are not trans- 
formed.) If the Fock space is endowed with its usual inner product (so that is adjoint to c\), the transformation 
is nonunitary in the bosonic sector since, if h,h^ are adjoints of each other, the transformed canonical pair b\ —6 no 
longer are. One can check that no mathematical problems are caused by this nonunitarity. After the transformation, 
the Fock space vacuum |0) and its dual (0| are annihilated by 

cValO) = cijO) = and (0|cV, - (0|cL, = 0. 

We may imagine that the "positive-energy states" {+a) are empty while the "negative-energy states" (—5) have been 
filled. Of course, for unitary bosons no such thing as a filled negative-energy sea exists, which is another way of seeing 
why the above canonical transformation must be nonunitary. Note also 
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? ^,Fci^,p|0} = |0) X Yl ^-/3.Bci^,BlO) = |0) X i-S^^n.). 

I3=l /3=1 

(The summation convention has been temporarily suspended for clarity.) 

Consider now the set of operators E^^g ■= c^a^b- commutation relations for the c, c these bilinears form a 

gl(nA^, nN) Lie superalgebra (n = n+ + 

Z^fci 1 p«i pfc' / ^\{\M + \B\){\C\ + \D\) rpkl T^ij 

Two subalgebras we shall have use for are the gl(A^) Lie algebra generated by the flavor-singlet operators J^a '^a'^A' 
and the gl(n,n) Lie superalgebra generated by the color-singlet operators J^i'^A'^B- These subalgebras commute 
and, moreover, they are maximal relative to each other, i.e. the maximal subalgebra of g\(nN,nN) whose elements 
commute with all elements of g^A'^), is gi(n, n), and vice versa. 

In the sequel, we focus attention on that part of Fock space where all gl(A^) generators vanish identically: 

c^c^ [flavor state) = 0. 

A 

This subspace contains the vacuum: 

E ^aC^IO) = E {^-0,B^-P,B + c-LftFci^,F) |0) = |0) X + n_) = 0, 

A 0=1 

and will be called the "flavor sector" (or color-neutral sector). We claim that the action of gl(n, n) on the flavor sector 
is irreducible^ i.e. every state of the flavor sector can be reached by a multiple action of the flavor operators c^c^ 
on the vacuum. The argument goes as follows. By employing the occupation number representation of Fock space, one 
easily sees that gl(nA^, nN) acts irreducibly on the subspace of Fock space selected by the condition ^ c^c^ = 0, 
which contains the color-neutral sector and might be called the "charge-neutral" sector. Thus every charge-neutral 
state \N) can be obtained by acting on the vacuum (or lowest- weight) state |0) with the raising operators, c!|.^c^j,, of 
gl(nA^, nN): 

If i? £ Gl(A^) is a rotation in color space, such a state transforms as a number of copies of the vector representation 
S'a ^ Y^j K^^c^A and the same number of copies of the co-vector representation i-^ Yji^~^y^'^A- Now assume \M) 
to be color-neutral, i.e. J2a ^a'^aI-^) ~ ^- order for this equation to hold, the indices of the expansion coefficients 
have to be contracted pairwise (tt here stands for permutations): 

Klb['::^X = E f{aibi...arbr;n)6,^^Q^y..6,^^Q^), 

by a standard result on the reduction of tensor-product representations of G1(7V). This form of the expansion 
coeflacients permits to express the operators that create |7V) out of the vacuum, in terms of the raising operators 
Tlii'^+a^^-b °f flavor superalgebra gl(n,n). Thus the vacuum |0} is a cyclic vector for the action of gl(n,ri) 
on the flavor sector. Because the action of every particle-hole creation operator ^ - c!|_^c!_f, can be undone by the 
corresponding annihilation operator c!_jc!^(j, the existence of a cyclic vector implies irreducibility. (q.e.d.) 

In addition to the Lie superalgebra gl(n,n) we shall need a corresponding Lie supergroup, G\{n\n). Recall that 
G :— Gl(n|n) is defined as the group of regular complex supermatrices of dimension (n + n) x [n + ?T.).n We associate 
with each element g G G an operator Tg on Fock space by 



^One advantage gained by the nonunitary canonical transformation c'_j, c^;,, c^j, ( — l)'''''''^c^i, is that it allows us to 
present all gl(nAi', nN) generators in the form cc, whereas in the original presentation they would have been be of various types 
cc, cc and cc. 

^The change in notation from the Lie superalgebra gl(n, n) (komma) to the Lie supergroup G\{n\n) (vertical bar) reminds us 
of the fact that the definition of the latter requires the specification of some parameter Grassmann algebra, whereas the former 
does not. 



15 



Tg = exp ^A(lng) 



AB'-B 

A,B,i 



In view of the multivaluedness of the logarithm, we need to demonstrate that this operator is well-defined. As the 
diagonalizable supermatrices are dense in G, it suffices to do so for an element 5 G G of the form g = SXS~^, where A is 
a diagonal matrix containing the eigenvalues of g. The eigenvalues of In (7 = ^(In X)S~^ are defined up to permutations 
and shifts by integer multiples of 27ri. We follow the convention that the Grassmann variables in the supermatrix 
Ing anticommute with the fermionic operators of the set {c^, c^}- It is then obvious that the transformation c\ |—^ 
7a = EB('S'~^)AsCi3, c^A.^ Ta = J2b^bSba is canonical. Note c^^iln g) ^gd^ = X^TAT^lnAA. Now the 
number operators = 7^7^ have integer eigenvalues, so the ambiguity InA^ ~> IuAa + 27rm is unobservable in 

Tg = exp (^j^ j 7a 7a ^^-^a^^ ^^"^ '^g is indeed well-defined. 

We next claim that the mapping g 1— > Tg is a homomorphism of Lie supergroups: 

TgTh = Tgh. 

The proof goes as follows. First, we use Toxpx to transform by 

Ca ^ T'expx Ca T^-^p^ = expad(4^i3cCc) ^a = CBlexpX)^^, 



where ad(A)_B := AB — BA stands for the commutator, and the summation convention is now back in force. On 
setting expX — g we get TgC\T~^ = (^bSba ^^^^ 

'^g'^h^A'^h ^ ~ '^g (^B^Ba) ^ " ^CdCB^^BA — ^Bi9h)BA — "^gh^ATgh , 

from which wc infer that the product T^f^TgTf^ commutes with c\. For identical reasons, T^j^TgTj^ commutes also 
with c^. Now, because the set of Fock operators {c^,c^} act irreducibly on Fock space, the vanishing commutator 
implies that T'f^TgTy^ is a multiple of the unit operator: 

T;HTgT^ = fi9,h)yil. 

The second step is to show }{g,h) = 1. For that we use a formula for the differential of the exponential map of 
any Lie algebra onto its Lie group: 



e^+'^ =e^Tx{X), where 

dt I t=o 



T- V^/ ad"(^) 1 



-&<i(X) 



(n + 1)! ad(X) 

This formula trivially extends to the case at hand, namely gl(n|n), a Lie algebra with Grassmann structure fs^ . The 
set where Tx has an inverse is dense in gl(n|n). We put g — expX, h{t) = expfF, and apply the formula for the 
differential twice, the first time in the reverse order: 



d 
'dt 



^^/cxp(x)cxp(ty)-^ 
~ 'dt 



^P^cxp(x+r-i(ty)) 
_^e^Vc\{X + T^\tY))^^c^B 

^ e^Vc\{T^oT^\tY))^^c^s, 



d 
di 



exp (c^^abCb) — 



T, 



d 



t=o 



4=0^^''^*^ dt 



TgThU). 



By integrating the last equation one arrives at the desired representation property, Tgh =^ TgTh. (q.e.d.) 
The subgroup H = Gl(n+|n+) x Gl(7i_|n_) C G consisting of elements h = diag(A, D) acts by 

7diag(A,D) = exp (c'_^„(ln A)„(,cVh + cLa(lnD)„f,cLfc) . 
This subgroup stabilizes the vacuum: 
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Tdiag(A,i3)|0) = |0)exp7V^(-l)l"l+i(lnD)aa - |0) SDetD- 



N 



With /i(diag(A, 13)) := SDetD ^ , we say that the vacuum carries a one-dimensional representation /i of H: 

Th\0) = \0)n{h), and {0\T-^ = ^i{hy^{0\ {h £ H). 

After all these preliminaries we consider the generalized coherent states Tg\0) [g G G) and, in particular, the 
operator P defined by the integral 

P = I DgHTg\0){0\T-\ 

(Recall that Dgn denotes the G-in variant Berezin measure of the coset space G/H, and the integration domain 
Mb X Mp was specified in Sec. [IB .) The integrand is a function on the coset space: 

TgH\0){0\T-,! = TgT^\0){0\T,;'T-' = rjO)(0|T-i {h e H), 
so the integral is well-defined. The following calculation shows that the operator P commutes with Tg^ for any go £ G: 

Tg„P = I DgH Tg„Tg\0){0\T^' 

= J Dg„ Tg„g\0){0\Tg-' 

= J DgHTg\0}{0\T-X„^PTgo, 



where the invariance of Dgn under translations g gog was used. (This invariance is crucial and is what determines 
D{Z, Z) = Dgn in (||), up to multiplication by a constant.) As was explained earlier, the action of the Lie superalgebra 
of G on the flavor sector is irreducible. By Schur's lemma then, the fact that P commutes with all operators Tgg leads 
to the conclusion that P is proportional to the identity on the flavor sector. 

Let us determine the constant of proportionality. To that end, let tt : G ^ G/H be the canonical projection 
assigning group elements g to cosets n{g) = gH. We make a decomposition g = s{T:{gyih(£), where s(7r(f/)) takes 
values in G and h{g) takes values in H. (Such a decomposition flxes a "choice of gauge" |2j,|l^, and s : G/H ^ G is 
called a (local) section of the bundle tt : G ^ G/H.) In this way we get 

T'glO) = Ts{-,^(g))Th(g)\^) = T,(^(3))|0) n{h{g)) , 
(0|r-i = mg))-' m;(l(g)y and T^|0)(0|r-i = T (.(,))|0)(0|T7i(^„. 

Following the main text we put 9 — ^) and parameterize cosets TT{g) by the pair {Z = BD^^^ Z = CA^^). For 
s we choose 



.(.(,))., s(Z,Z)^( ^]-^Xu. T'~^Tui' 

^ ^ ' ' \ z(\- zzy^i^ (i-zz)-^/^ 

1 Z\ /(l-ZZ)+i/2 V ^ 

1 j 1^ (1-ZZ)-V2 ) [z 1 

The first line shows that this is a valid choice of section, i.e. we indeed have Z = BD^^ and Z = CA^^. Using the 
second line to translate s{Z, Z) into an operator on Fock space, we obtain 



|0) ^ cxp {c\^Z^^c'_^) exp (^ic^^ ln(l - ZZ)^^c\f, 

-i? ,ln(l - ZZ),,cL,) exp (cL,Z,,cV,) |0) 
= exp {c\,Z,,c'_,) |0) SDet(l - ZZ f/' \Z), and 
^^Ki,z) = (0|^.(-z,-z) - SDet(l - ZZ)^/2(0|exp {^ctJa^c^^ {Z\. 
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The expression for the invariant measure Dgn in the coordinates Z, Z is D{Z, Z), so 

f Dgn r,|0)(0|T-i = / D{Z,Z) T ^Mm;' = f D{Z, Z)\Z){Z\. 



' Mb X J\/f 

Taking the vacuum expectation value of P we find 



(0|P|0) = J D{Z, Z) {Q\Z){Z\Q) = j D{Z, Z) SDet(l - ZZ)^ = j Dhn{Z, Z) = 1, 



by our choice of normalization, see Sec. II A . It follows that the constant of proportionality between P and the unit 
operator on the flavor sector is unity. Because {Z\ vanishes on states that are not color-singlets, P annihilates all 
states outside the flavor sector. Hence P projects Fock space onto the flavor sector. 
The final ingredient we need are the Bose-Fermi coherent states 

exp (c;>;, + v^LfcC^b) |o). 

Note that the ordering in the exponential matters as the Grassmann variables are taken to anticommute with the 
fermionic operators: (?^g^ip+a — (— l)^°^'i/'lnc!|_„, as before. These Bose-Fermi coherent states span the entire Fock 
space (or, rather, the "Grassmann envelope" [ p8[ thereof). They can be projected on the flavor (or color- neutral) 
sector by making a unitary rotation in color space, c\ i-^ c'^U^^, i-^ U^^cf^, and averaging over all such rotations. 
Therefore the operator P that projects on the flavor sector acts on Bose-Fermi coherent states as 

Pexp (cV^V^V, + Cfccl,) |0) = / dU exp (c7^V+„V^V, + ^Itc'^.W') |0). 
With all these tools in hand, the statement (0) is proved by the following computation: 



D^i^{Z, Z) exp (^V.Z^.Cfc + i>UZbJAa) 
D^^N{Z, Z) (0| exp (V^^.cV, - c^^U) exp (cVa^afec'l,) |0) 

X (0| exp (-cL^Z^^d^^ exp (cVaV'^a + i^-b<^b) |0) 
= (0|exp(V;V,cV<, bCb) (/ D{Z,Z)\Z){Z\^ e^v {^+a4'\a + i^U^-b) |0) 
= (0| exp (v:;V,c;, - c^.C^) ^exp (c;,V;a + ^U^-b) |0) 

dU (0| exp (^V,cV„ - cLbV-.b) exp U+JJ^'i^^, + ^L.t/^Vl,) |0) 

du exp (^;,c/'Vi, + ^L,C7^'>i, 



V{N) 



V(N) 

The first equality sign is an elementary consequence of the commutation relations of the Fock operators c, c and their 
action on the vacuum |0). The second and third equality signs recognize the projector P on the flavor sector. The 
fourth equality sign implements P by averaging over all unitary rotations in color space. The last equality sign is 
elementary again. 
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